Electronic correlation effects and the Coulomb gap at finite temperature 
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We have investigated the effect of the long-range Coulomb interaction on the one-particle excitation 
spectrum of n-type Germanium, using tunneling spectroscopy on mechanically controllable break 
junctions. The tunnel conductance was measured as a function of energy and temperature. At low 
temperatures, the spectra reveal a minimum at zero bias voltage due to the Coulomb gap. In the 
temperature range above 1 K the Coulomb gap is filled by thermal excitations. This behavior is 
reflected in the temperature dependence of the variable-range hopping resitivity measured on the 
same samples: Up to a few degrees Kelvin the Efros-Shkovskii InJ? tx T -1 / 2 law is obeyed, whereas 
at higher temperatures deviations from this law are observed, indicating a cross-over to Mott's 
ln_R oc J 1 " 1 / 4 law. The mechanism of this cross-over is different from that considered previously in 
the literature. 

PACS numbers: 74.50. +r, 74.40. +k, 74.60. Jg, 73.40. Jn 



The electronic density of states (DOS) near the Fermi 
level is an important physical quantity for understanding 
electrical transport mechanisms in strongly localized sys- 
tems like impurity bands in doped semiconductors. 
They consist of sites with random positions and random 
energies . The electron wavefunction is localized with 
a localization radius a smaller than the average nearest- 
neighbor distance between sites which is of order A -1 / 3 
(A = concentration of impurities). At low temperatures 
the electrical resistivity of such systems is governed by 
variable-range hopping (VRH) , which means that the ac- 
tivation energy for a hopping process decreases continu- 
ously with temperature 0||. The role of the inter-site 
Coulomb interaction between electrons in the hopping 
regime was first addressed by Pollak and by Srinivasan 
H . They showed that in localized systems the Coulomb 
interaction creates a deep depletion of the one-particle 
DOS near the Fermi energy Ep. Efros and Shklovskii 
(ES) called this depletion 'Coulomb gap' and showed that 
the DOS near E-p varies as g(E) — go (E — Ep) _1 for di- 
mensionalities D — 2 and D — 3, respectively |^||. This 
leads to a VRH hopping resistance R oc exp (Tpg/T) 1 ^ 2 
§,0, in contrast to Mott's R oc exp (T M /T) 1/4 law § for 
which a constant DOS g(E) — go is assumed. 

A cross-over between these two temperature laws has 
been expected theoretically H and also observed exper- 
imentally ||[|. According to the traditional interpreta- 
tion of this cross-over the energy range of the phonon- 
assisted tunneling (hopping) becomes larger than the 
width of the Coulomb Gap Aqg above the cross-over 
temperature. In this case the Coulomb gap does not af- 
fect the hopping resistance, thus resulting in Mott's law. 
However, another kind of crossover from an ES-type to 
Mott-type temperature law may occur, namely the fill- 
ing of the Coulomb gap by thermal excitations. Such a 



filling of the Coulomb gap with increasing temperature 
was predicted by Monte-Carlo simulations of a classical 
Coulomb glass []10| pL2[1 . More recently, it was observed by 
tunneling sprectroscopy on Si:B samples near the metal- 
nonmetal transition. The temperature dependence of the 
DOS should, in turn, affect the temperature dependence 
of the bulk resistivity. Surprisingly, except for Ref. [fl2f , 
this effect has not yet been taken into account in the 
literature. 

In this letter we present evidence that in insulating 
doped Ge the second mechanism is verified. We show 
by comparing the temperature dependence of the DOS, 
derived by tunneling spectroscopy, with that of the VRH 
resistivity that the observed deviations from the ES 
law above 1 K correspond to the thermal filling of the 
Coulomb gap and not to the traditional cross-over mech- 
anism. 

The electronic DOS in solids can be dicectly probed by 
tunneling spectroscopy and photoelectron spectroscopy, 
for example. But the small size of the width of the 
Coulomb gap Aqg strongly restricts the useful spectro- 
scopic techniques. At present only the tunneling spec- 
troscopy has the required energy resolution. Massey and 
Lee were the first to directly observe the Coulomb gap 
in a doped uncompensated semiconductor Si:B [pL 3| pL 5|| - 
They used planar tunnel junctions between boron-doped 
Si samples and a Pb counter-electrode with an insulating 
dielectric as barrier. The superconducting quasi-particlc 
DOS of the lead electrode has been observed, proving 
quantum tunneling. Suppressing superconductivity of 
lead in a magnetic field of B = 200 mT allowed then 
to measure the DOS of the Si:B electrode against the 
constant DOS of normal lead. However, planar tunnel 
junctions are difficult to prepare, especially with germa- 
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As alternative method we have proposed recently that 
it is possible to realize tunneling across a semiconductor 
break junction due to the lateral confinement of small 
point contacts |L6]]. This is a well established technique 
to investigate superconductors and metals, both in the 
regime of direct metallic contact and tunneling across 
a vacuum barrier |l7|. To ensure a voltage drop con- 
fined to the junction itself the contact resistance is set 
to a value larger than about 10kf2 at 1 K. Our junctions 
have lateral dimensions of less than about 100 nm. The 
hopping length amout to about 150 nm at 1 K increas- 
ing towards lower temperatures. Charge is transported 
across the junction by a single hopping event at a rather 
well-defined energy. This justifies our calling the trans- 
port tunneling and to interpret the voltage drop as an 
excitation energy. On the other hand, the contact diam- 
eter must be large enough to inhibit the formation of a 
depletion layer. Such a layer adds an additional tunnel- 
ing barrier |^8| , and it may also vary the local DOS in 
the contact region. We estimate a ~ 10 nm lower bound 
for useful junctions. 

Our undoped samples were grown by the Chochral- 
skii method using highly enriched (up to 93 % 74 Ge) 
germanium. Neutron-transmutation doping (NTD) en- 
sured excellent homogeneity of 75 As as shallow donors 
Pj . The main nuclear reaction of Ge with thermal neu- 
trons is 74 Ge (n, 7) 75 Ge — ► 75 As. As by-product, a small 
fraction of 75 Ga is produced as acceptors. This gives a 
compensation of K — Nc, a /NAs — 12 %, and sets the 
Fermi level inside the impurity band. The donor concen- 
tration Nd = Nas was below, but close to, the disorder- 
driven metal-insulator transition of Ge at a critical im- 
purity concentration of 3.4 • 10 17 cm~ 3 p9| . 

For our tunneling experiments the samples were cut 
into 1 x 1 x 10 mm 3 slabs with a 0.5 mm deep groove to 
define the break position within the (111) cleaving plane 
of germanium. The samples were glued onto a flexible 
bending beam, electrically isolated but thermally well 
coupled to the cold plate. They were broken at low tem- 
peratures in the ultra-high vacuum chamber of the refrig- 
erator. The contact size could be adjusted in situ with a 
micrometer screw and a piezo tube. For further details of 
the setup see Ref. po| . The dl/dU spectra of junctions 
with small resistance (less than about 100 kfi) were ob- 
tained by means of the standard four-terminal method 
with current biasing. The current-voltage characteris- 
tics of junctions with high resistance (larger than about 
100 kf2) was recorded using the standard two-terminal 
method with voltage biasing. In the latter case the bulk 
samples contributed at most 5 % to the total resistance. 

All investigated Ge break junctions have rather similar 
characteristics. Fig. |l| shows typical spectra of a sam- 
ple with Nd = 1.26 x 10~ 17 cm~ 3 as function of voltage 
and temperature. The spectra have a pronounced mini- 
mum at low temperatures. We believe that this anomaly 
represents the Coulomb gap. Between 100 mK and 1 K 



the spectra depend only weakly on temperature. Above 
about T = 1K the Coulomb gap becomes filled by ther- 
mal excitations. It has almost vanished at T = 6K. 

In order to investigate how this temperature depen- 
dence of the DOS affects the hopping resistivity we have 
measured the resistance of the bulk sample using the 
standard four-terminal technique. Fig. shows the re- 
sistance as function of T -1 / 2 and T^ 1 ™, respectively. 
At low temperatures lni? oc T -1 / 2 , as expected for the 
ES law. The resistance deviates from this behavior at 
T -1 / 2 < 0.5 K- 1 / 2 or at T > 4K, nearly coincidig with 
the temperature at which the Coulomb gap is suppressed 
according to the tunnel data, see Fig. [|. 

We have also considered the traditional scenario of the 
crossover between the ES and Mott law in which one as- 
sumes that the hopping energy range becomes broader 
than Coulomb gap || . ^From an analysis similar to that 
of EJ for the traditional cross-over mechanism we ob- 
tained a width the Coulomb gap of Aqg = 0.25 meV 
which is considerably smaller than the observed Aqg = 
2.5 meV. But the latter value would correspond to a much 
higher crossover temperature than observed experimen- 
tally. 

We turn now to a more detailed discussion of the tem- 
perature dependence of the DOS. To extract the DOS 
from the spectra we first tried to remove the energy- 
dependent part at high voltages. Several possibilities 
were tried, with only slight variation of the final re- 
sult. According to Ref. [ p^| , this high-energy tail can 
be roug hly described by g{E) oc 1 + y/{E - E F )/5. In 
this model, the parameter S represents a correlation en- 
ergy, which is almost independent of impurity concentra- 
tion. For our experiments 6 ~ lOmeV, a rather large 
value when compared to the results for Si:B |^5|. Al- 
ternatively, a Schottky-type behaviour was used with 
dl/dU oc exp (U/Uqq). The parameter Uqo may then 
represent the properties of an additional barrier due to 
the depletion layer. As we do not know which of the 
two possibilities is correct, we normalize the spectra at 
low temperatures with respect to that corresponding to 
the highest temperature. The shape of those normalized 
spectra is almost flat outside the Coulomb-gap anomaly. 
They are then fitted using 



0(£,T)oc7(T)+[l- 7 (T)] 



\E-Eg 



(Acg{T)/2)* + \E-E f \ 



(1) 



The parameter 7 describes a 'residual' DOS at the Fermi 
level and Aqg is the width of the Coulomb gap (FWHM) . 
The DOS derived by ES is recovered when 7 = and 
s = 2g|. 

The experimental DOS of our samples strongly devi- 

and 
IL3ML5I. 



ates from the simple square-law derived by ES [p 
which was also found experimentally for Si:B 
Taking into account the Fermi distribution and the DOS 
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Eq. |l| on both sides of the junction, our analysis yields 
for the sample far from the metal-nonmetal transition 
s = 3. Fig. H shows how Acg and 7 depend on tem- 
perature. Both saturate at low temperatures. Analytical 
as well as numerical simulations for nonmetallic disor- 
dered systems have predicted several different relation- 



ships: power laws g(E, T = 0) x(E-E F §J2| and 



g(E, T = 0) oc {E-E F ) 2 - 7±0A |lfj as well as an exponen- 
tial dependence g(E,T = 0) oc exp[-(A/(E - £ F )) 1/2 ] 
p4j . A power law with s = 3 is consistent with the results 
of Ref. 0. 

The temperature dependence of the zero-bias DOS also 
reveals a power law N(Ep,T) oc T x with an exponent 
x = 0.8. This differs from x = 2.7 derived by |Tc|] , but it 
agrees quite well with x = 1 obtained by the simulations 
of Ref. 0. 

Because of the experimental DOS, especially its strong 
temperature dependence, the observed ES-type be- 
haviour of the bulk resistivity needs an explanation. We 
propose the following scenario to interpret our results: i) 
At low temperatures the Coulomb gap ist temperature 
independent. This corresponds to the ES regime of the 
VRH conductivity; ii) the power law behaviour of the 
residual DOS j(T) causes the transition from the ES to 
the Mott law; Hi) at high temperature the DOS is con- 
stant and one expects a Mott-type VRH conductivity. 
The mechanism of the transition between ES- and Mott- 
type behaviour is therefore completely different from that 
considered in Refs. p|,^H 

To summarize, the tunnel conductance of small break 
junctions of our germanium samples shows a minimum 
of the DOS near the Fermi level. This minimum repre- 
sents the Coulomb correlation gap. Up to about 1 K, the 
width of this anomaly depends only weakly on T. This 
is the ES-regime. Above about 1 K, the anomaly smears 
out. Consequently deviations from the ES law occur. 
Above about 6K, the anomaly has vanished, and the 
DOS becomes constant. The sample is then in the Mott 
regime. According to our interpretation, the cross-over 
between ES and Mott law at such low temperatures is not 
due to the temperature-enhanced range of hopping ener- 
gies. But it origins from the suppression of the Coulomb 
gap by thermal excitations. These break-junction tun- 
neling results are consistent with the electrical resistivity 
of the same samples, which also shows a rather wide tem- 
perature range for the transition between ES and Mott 
variable-range hopping. 

Acknowledgments: This work was supported by the 
German SFB 252 Darmstadt/Frankfurt/Mainz and the 
Russian RFFI 97-02-18280. 



[1] P. W. Anderson, Phys. Rev. 109, 1492 (1958). 

[2] N. F. Mott and E. A. Davis: Electron Processes in Non- 
crystalline Materials, 2nd ed. (Clarendon, Oxford 1979) . 

[3] B. I. Shklovskii and A. L. Efros, Electronic properties 
of doped semiconductors, (Springer- Verlag, Berlin, New 
York 1984). 

[4] M. Pollak, Disc. Faraday Soc. 50, 13 (1970); M. Pollak, 

Proc. Roy. Soc. (London) A325, 383 (1971); M. Pollak, 

Phil. Mag. 23, 519 (1971). 
[5] G. Srinivasan, Phys. Rev. B 4, 2581 (1971). 
[6] A. L. Efros and B. I. Shklovskii, J. Phys. C 8, L49 (1975). 
[7] A. L. Efros, N. V. Lien, and B. I. Shklovskii, J. Phys. C 

12, 1869 (1979). 
[8] A. N. Ionov, I. S. Shlimak, and M. N. Matveev, 

Solid State Commun. 17, 314 (1983). 
[9] I. S. Shlimak, L. J. Zaribin, A. N. Ionov, V. M. Vorobkalo, 

A. G. Zabrodskii, and I. Yu. Nemish, JETP .... [Pis'ma 

v JTPh 1983, 9, 877 (1983)] 
[10] M. Sarvestani, M. Schreiber, and T. Vojta, Phys. Rev. B 

52, R3820 (1995). 
[11] E. R. Grannan and C. C. Yu, Phys. Rev. Lett. 71, 3335 

(1993). 

[12] I. Shlimak, M. Kaveh, R. Ussyshkin, V. Ginod- 
man, S. D. Baranovskii, P. Thomas, H. Vaupel, and 
R. W. van der Heijden in: Hopping and Related Phe- 
nomena, O. Millo, Z. Ovadyahu, Eds., p. 50 Jerusalem, 
1995). 

[13] I. G. Massey and M. Lee, Phys. Rev. Lett. 75, 4266 

(1995) . 

[14] J. G. Massey and Mark Lee, Phys. Rev. Lett. 77, 3399 

(1996) . 

[15] M. Lee, J. G. Massey, V. L. Nguyen, and B. I. Shklovskii, 

Phys. Rev. B 60, 1582 (1999). 
[16] B. Sandow, K. Gloos, Yu. Naidyuk, R. Rentzsch, and 

A. N. Ionov, phys. stat. sol. (b) 205, 281 (1998). 

[17] C. J. Muller, J. M. van Ruitenbeek, and L. J. de Jongh, 

Physica C 191, 485 (1992). 
[18] B. Sandow, K. Gloos, R. Rentzsch, and A. N. Ionov, 

Physica B, in print. 
[19] A. N. Ionov, M. J. Lea, and R. Rentzsch, JETP Lett. 54, 

470 (1991) [Pis'ma Zh. Eksp. Teor. Fiz. 54, 470 (1991)]. 
[20] K. Gloos and F. Anders, J. Low Temp. Phys. 116, 21 

(1999). 

[21] I. S. Shlimak, M. Kaveh, M. Lea, P. Fozooni, P. Stefanyi, 

and A. N. Ionov, Physica A 200, 483 (1993). 
[22] S. D. Baranovskii, A. L. Efros, B. L. Ge'mont, 

B. I. Shklovskii, Solid State Commun. 27, 1 (1978). 
[23] S. D. Baranovskii, A. L. Efros, B. L. Ge'mont, 

B. I. Shklovskii, Solid State Commun. 27, 1 (1978). 
[24] A. L. Efros, J. Phys. C 9, 2021 (1976). 



and 



and 



FIG. 1. dl/dU vs. U spectra of the Ge sample with 
N d = 1.26 ■ 10 17 cm -3 . 



FIG. 2. Electrical resistance of the bulk Ge sample R 
vs. T 1/4 and T 1/2 , respectively. N d = 1.26 ■ 10 17 cm -3 . 
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FIG. 3. Coulomb gap Acg and residual DOS 7 vs. tem- 
perature T. N d = 1.26 • 10 17 cm~ 3 . 
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